The proceedings of the Physics of Accreting Compact Binaries 



1 



CYCLOPS-X: Simultaneous optical and X-ray modeling of 
polars 



Karleyne M. G. Silva 1 , Claudia V. Rodrigues 1 and Joaquim E. R. 
Costa 1 

(1) Divisao de Astrofisica, Instituto Nacional de Pesquisas Espaciais, Av. dos 
Astronautas 1758, 12227-010, Sao Jose dos Campos-SP Brazil 



Abstract 

From the optical to the X-ray frequencies, most of the continuum emission 
from AM Her systems originates in the post-shock region. Hence, using multi- 
spectral data can be useful to restrict the physical and geometrical properties of 
that region. In spite of that, distinct codes are used to model these frequency 
ranges. CYCLOPS is a code to model the optical cyclotron emission of polars. In 
this contribution, we present the first version of CYCLOPS-x, an improvement of 
CYCLOPS to fit simultaneously optical and X-ray data. The new code adds the 
bremsstrahlung as an emission process as well as the effects of X-ray absorption by 
the upper portion of the accretion column. As a first application of CYCLOPS-X, 
we present X-ray light curves using two sets of parameters provided by the optical 
modeling of CP Tuc. These two cases have very similar optical emission, but 
quite distinct light-curves in high energies. It illustrates the need of simultaneous 
modeling of the optical and X-ray emission to a proper description of the magnetic 
accretion in polars. 

1. Introduction 

Polars are a subclass of the cataclysmic variables (CVs) in which the accreted 
material is channeled by the magnetic field of the white dwarf (WD) forming an 
accretion column. Near the WD surface the material forms a shock front increasing 
temperature and densities. The bulk of the emission comes from the base of the 
column also called post-shock region. The mean emission processes are cyclotron 
(optical and infrared) and bremsstrahlung (X-ray). The radiative cooling results 
in a structured region with temperature and density varying as a function of the 
height of the post-shock region. 

The X-ray emission of polars is commonly reproduced combining shock models 
and plasma emission models, as Mcka model [T] available in the xspec Spectral 
Fitting Package [2]. First, the post-shock region is divided in layers and each 
layer has a temperature and density given by the post-shock structure. Then, the 
cmissivity is calculated for each layer and the total emission is the sum of the 
plasma cmissivity of all layers, assuming an optically thin region [3]. 

The X-ray fittings are usually done without considering the optical data, in 
part because distinct codes are used to model each spectral range. An exception 
is the modeling of polar RX J2115-5840 [4]. The cyclotron 2D emitting region, 
obtained from the fit of optical data, is used to define the location of the X-ray 
emission. However, the density and temperature profiles adopted in fitting each 
energy band are not the same. 

The absorption due the pre-shock flow located above the post-shock region is 
important in high energies. Customary, the absorption, which varies with orbital 
phase, is arbitrarily chooscn in order to fit the observational data. 

The CYCLOPS code solves the radiative transfer of the cyclotron process in the 
post-shock region adopting a 3D treatment and specific functions to describe the 
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variation of temperature and density due the radiative cooling. The code provides 
the four Stokes parameters which can be used to reproduce optical light and 
polarimetric curves from polars [5] . Here we present the modification of CYCLOPS 
to include bremsstrahlung emission and to prepare the code to simultaneously 
model the optical and X-ray continuum emission of post-shock regions of polars. 

2. CYCLOPS-X 

CYCLOPS-X has exactly the same input parameters of CYCLOPS [5]. Hence, 
no additional parameters are needed to model X-ray data. The post-shock and 
pre-shock regions are defined by the magnetic field lines and divided in several 
elements of volume according to the chosen spatial resolution. The emergent flux 
in a given orbital phase is the sum of the fluxes coming from all lines of sight. 
The code solves the radiative transfer to each line of sight in steps, considering an 
inhomogcneous post-shock region. This naturally incorporates the shock structure 
and optical depth effects. 

The emissivity of the bremsstrahlung is calculated in each voxel according to 



is the temperature and g is the non-relativistic Gaunt factor [I] . 

The absorption is considered only in the phases in which the pre-shock region is 
between the observer and the post-shock region. This is done for each line of sight 
that composes the total flux in a given orbital phase. The extinction processes 
considered are photoelectric absorption [7] (a p h,) and Thomson scattering (a t h)- 
The flux observed in each orbital phase (F p h a ) is sum of the emergent flux in a 
line of sight (Fq) attenuated by the upper part of the column according to: 



where s is the optical path in the line of sight and / is the number of lines of sight 
used to represent the region. 

3. An example of CYCLOPS X-ray light curves: CP Tuc 

In order to illustrate the light curves produced by CYCLOPS-X, we chose the 
polar CP Tuc. 

Initially, it was suggested that the post-shock region of CP Tuc is always visible 
and that the modulation observed in the X-ray light curves would be due absorp- 
tion by the upper accretion column in specific orbital phases [5]. After that, the 
modeling of the optical light and polarization curves indicated a different geome- 
try to CP Tuc, in which its post-shock region is self-eclipsed (by the WD) during 
half of the orbital period generating the modulation [9] . Recently, our group pre- 
sented new data in two optical bands. CYCLOPS modeling of these data shows 
that it is possible to find good models in the two cases above: the absorption and 
the self-eclipse scenarios. (Please see Cyclotron modeling of the polar CP Tuc, 
Rodrigues et al., this meeting, for more details.) 

We have used CYCLOPS-x to obtain the X-ray light curves in each case using 
the parameters found by the optical fits. In doing that, we are considering that 
optical and X-ray data are representative of a same brightness state. Figure [TT1 
presents the observed [8] and the modeled light curves in three bands: 0.7-2.3 keV 
(top), 2.3-6.0 keV (middle) and 6.0-10.0 keV (bottom). 
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Fig. 1. Comparison between X-ray light curves of CP Tuc (points) and the modeled 
X-ray light curves. The abcissa is orbital phase and runs from to 2. Left: Ab- 
sorption scenario (solid line= model with absorption; dashed line= model without 
absorption). Right: Self-eclipse scenario. 



Absorption scenario: Figure shows 11 images in different orbital phases 
produced by CYCLOPS-x. Specifically, it is shown the X-ray images of the post- 
shock region with (top) and without (bottom) the pre-shock absorption. Figure [TT1 
shows that the absorption reproduces the energy dependence of the modulation, 
but the amplitude was insufficient to fit the data in the hardest bands. As the 
density of the pre-shock region has an upper limit defined by the minimum density 
of the post-shock region, an improvement of the X-ray fit may be achieved by an 
increase of the density in the post-shock region. However, it is necessary to check 
the influence on this on the optical modeling. 

Self-eclipse scenario: Figure [3] presents 11 images of the emitting region 
distributed along the orbital period. In this case, the absorption of the pre-shock 
region is not important, because it never is the line of sight. This model does 
not reproduce the dependence with energy as showed in Figure [TTJ How could 
this be improved? In this specific model, only the base of the accretion column is 
self-eclipsed. That region is responsible by most of the emission in soft energies, 
according to Equation [1] and the post-shock structure. A possibility to improve 
the X-ray fit is to change the temperature and the height of the post-shock region. 
A higher column could present a larger range of temperatures and densities: it 
could provide a modulation more pronounced in the soft band than in the hard 
band. 

In both cases, the sets of parameters that produce good models in optical bands 
do not reproduce the X-ray data. Nevertheless some important features of the data 
could be reproduced. 




Fig. 2. Images of the post-shock region in the absorption scenario. Up: considering 
absorption. Bottom: with no absorption. The images have different lookup tables. 





Fig. 3. Images of the post-shock region in the self-eclipse scenario. 



4. Conclusions 

We present the CYCLOPS-x, a code that calculates the optical and X-ray emis- 
sion of post-shock regions of polars. This code includes the absorption caused 
by the pre-shock region in a manner consistent with the geometry defined by the 
magnetic lines. We used the CYCLOPS-X to calculate the X-ray light curves of 
the two good optical models for CP Tuc. The optical emission of the two models 
are very similar, what is not true for the X-ray fluxes. This simple application 
shows the importance of consider as many data as possible to understand the 
post-shock region of a polar. We are planing to incorporate in CYCLOPS-x the 
simultaneous fit of X-ray and optical data and, then, apply this tool to the CP 
Tuc multiwavelength data. 
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